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LONG-TERM  GOALS 

Predictions  of  the  ocean  dynamics,  sediment  transport,  pollutant  dispersal  and  biological  processes 
require  knowledge  on  the  characteristics  of  turbulence  in  the  bottom  boundary  layer.  Modeling  of  the 
turbulence  requires  data  for  development/validation  of  closure  models.  Consequently  our  goal  is  to: 

a.  Measure  the  Reynolds  stresses  (free  of  wave  contamination),  velocity  profile,  dissipation  rate,  and 
turbulent  spectra  in  the  coastal  bottom  boundary  layer  using  Particle  Image  Velocimetry  (PIV). 

b.  Quantify  the  temporal  variation  of  turbulent  stresses,  production,  dissipation  and  buoyancy  flux  in 
relation  to  the  oceanographic  parameters  that  represent  the  local  environment,  such  as  waves,  currents, 
stratification,  internal  waves  and  the  nature  of  the  water-sediment  interface.  The  conclusions  will 
quantify  the  relative  importance  of  different  mechanisms  that  control  the  flow  and  turbulence. 

c.  Study  the  mechanisms  and  extent  of  sediment  re-suspension  process  by  simultaneously  measuring 
the  flow  structure  and  particle  distributions. 

d.  Use  oceanic  PIV  data  for  addressing  Sub-Grid  Scale  Modeling  issues  for  Large  Eddy  Simulation  in 
oceanic  flows. 

e.  Examine  the  structure  of  the  flow,  vertical  vorticity  transport,  formation  and  upward  migrations  of 
large  coherent  vortex  structures.  Presently  there  is  very  little  information  on  the  dynamics  and  impact 
of  large  coherent  structure  in  the  bottom  boundary  layer  on  turbulence  and  sediment  entrainment. 

OBJECTIVES 

Our  effort  for  the  past  year  focused  on  several  objectives: 

a.  A  field  trip  (funded  in  part  by  the  present  ONR  grant  and  in  part  by  NSF)  to  the  vicinity  of  LEO- 15 
took  place  in  the  fall  of  2001 .  We  have  used  a  new  platform  that  has  a  profiling  range  that  extends 
from  very  close  to  the  bottom  up  to  10  m  above  the  bottom. 
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b.  Analysis  of  the  PIV  data  obtained  in  several  deployment,  including  data  recorded  in  the  near  LEO- 
15  during  summers  2000  and  Fall  2001,  and  in  the  mouth  of  the  Delaware  Bay  in  2000.  In  these 
experiments  we  recorded  data  in  two  sample  areas  simultaneously.  The  data  has  been  used  for: 

1 .  Measuring  the  mean  velocity  distributions,  RMS  values  of  velocity  fluctuations  (free  of  wave 
contamination  by  using  structure  functions),  turbulence  energy  spectra,  dissipation  rates,  and  the 
Taylor  micro-scale  Reynolds  number.  The  results  are  used  for  characterizing  the  turbulence  in  the 
vicinity  of  LEO- 15  under  various  conditions  of 
mean  current  and  wave  amplitude. 

2.  Examining  the  flow  structure,  including  the 
dynamics  and  distributions  of  large-scale  vortices. 

3.  Measuring  the  Sub  grid  Scale  (SGS)  stresses 
needed  for  applications  of  Large  eddy  Simulations 
(LES)  for  modeling  of  oceanic  turbulence.  The 
results  are  then  used  for  estimating  the  SGS  energy 
flux  (dissipation)  and  comparing  it  to  the  viscous 
dissipation  rate.  The  SGS  dissipation  is  a 
quantitative  measure  of  energy  transfer  between 
turbulence  at  different  length  scales.  Thus,  its 

dynamics  sheds  light  on  the  interaction  between  Figure  1.  The  oceanic  PIV  system. 

structures  of  different  sizes. 


4.  Evaluating  the  performance  of  several,  commonly 
used  SGS  stress  models,  including  the  Samgorinsky, 
Eddy  Viscosity  model,  and  the  non-liner  model. 

APPROACH 

The  oceanic  PIV  system  is  described  in  refs.  12-13. 
The  present  system  has  improved  substantially  from 
the  original  setup  described  in  refs.  2  &  4.  A 
schematic  of  the  submerged  components  is  shown  in 
Figure  1.  The  laser  (a  pulsed  dye  laser  generating 
pairs  of  2  ps  pulses)  is  located  on  the  ship  and  the 
light  is  transmitted  through  optical  fibers  to 
submerged  probes.  Images  are  acquired  using  two 
2k  x  2k  CCD  cameras,  each  with  a  maximum 
sampling  rate  of  4  frames/s.  Each  camera  and 
associated  light  sheet  can  be  aligned  independently, 
in  the  same  or  different  planes,  near  each  other  or 
apart.  The  data  is  recorded  on  mass  data  acquisition 
systems  (one  for  each  camera),  allowing  continuous 
sampling  at  0.5  Hz  for  over  15Hrs.  The  data  analysis 
procedures  are  described  in  refs.  12,  13,  17  and  18. 


U  =  12.4  cm/s,  W  =  -0.2  cm/s 


-15  -10  -5  0  5  10  15 
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Figure  2.  Velocity  vector  map.  The  area 
mean  velocity  is  subtracted  from  each 
vector.  The  vertical  coordinates  represent 
the  actual  distance  from  the  bottom. 
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Each  instantaneous  velocity  distribution  consists  of  63x63  vectors,  and  the  vector  spacing,  8,  is 
indicated  in  Table  1.  A  sample  is  presented  in  Figure  2. 

The  submersible  components  of  the  PIV  system  are  mounted  on  a  seabed  platform,  which  can  be 
rotated  to  align  the  sample  areas  with  the  mean  flow  direction,  and  extend  vertically  to  sample  at 
different  elevations.  The  platform  is  a  5-stage  telescopic  hydraulic  cylinder,  with  a  vertical  range  of 
9.75m.  The  system  also  contains  a  CTD,  transmissometer,  dissolved  oxygen  sensor,  precision  pressure 
transducer,  clinometer,  digital  compass  and  video-microscope  for  sampling  the  particle  distributions  at 
high  magnification. 

WORK  COMPLETED 

In  refs  12  &  13  we  select  and  compare  several  data  sets  obtained  at  different  depth  to  represent 
conditions  of  relatively  high,  intermediate  and  weak  mean  flows.  They  also  represent  substantially 
different  turbulent  Reynolds  numbers.  Table  1  summarizes  the  selected  data  series.  To  characterize  the 
mean  flow  and  amplitude  of  surface  waves,  we  average  the  velocity  over  the  vector  map  to  obtain  the 
instantaneous  average  velocity  components  (U  and  W).  The  mean  current  is  characterized  by  U  and 
W  ,  the  overall  average  (over  all  distributions),  whereas  Ums  and  Wms  are  their  RMS  values, 
representing  mostly  the  effect  of  surface  waves  but  also  turbulence  at  scales  larger  than  the 
instantaneous  distributions.  The  data  of  Runs  A  and  B  were  obtained  at  the  mouth  of  the  Delaware 
Bay,  at  a  mean  near  the  bottom  velocity  of  38  cm/s  with  little  wave  motion.  Runs  C  to  F  took  place  in 
Fall  2001  near  FEO-15.  The  current  in  this  region  is  moderate  to  weak,  but  it  is  exposed  to  oceanic 
swell.  For  Runs  C  and  D  the  amplitude  of  the  wave  induced  velocity  is  of  the  same  magnitude  as  the 
mean  flow.  For  Runs  E  and  F,  the  mean  current  is  very  low,  namely  the  flow  consists  almost  entirely 
of  wave  induced  motion.  Data  were  collected  continuously  for  periods  of  20  min,  at  2  -  3.3  Hz,  and  at 
elevations  up  to  8.5  m  above  the  seabed.  Here  we  use  data  obtained  at  mean  elevations  (vector  map 
center)  of  0.55  and  about  2.5  m  above  the  bottom. 

Table  1.  The  selected  data  sets  showing  mean  velocity,  wave  induced  motion  (Ums,  Wms),  RMS  values 
of  turbulent  velocity  fluctuations  (u  ’.  w  ’),  Taylor  microscale  and  resulting  Re‘‘ . 
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Figure  2  demonstrates  the  presence  of 
large-scale  vortical  structures,  which 
appear  intermittently  in  all  the  data  sets 
except  for  those  with  very  low  mean 
velocity  (E  and  F).  Figure  3  presents 
mean,  one-dimensional  energy  spectra  of 
u  (En)  and  w  (E:t)  integrated  along  the 
streamwise  (k,)  directions.  These  spectra 
do  not  involve  use  of  the  Taylor’s 
Hypothesis.  In  addition,  each  graph 
contains  an  insert  containing  the 
distributions  of  Erf'kf'  En(k^  and 
(3  /  4)szjf 3 kf 3 E..  (kx ) .  As  is  evident,  the 
turbulence  is  anisotropic  at  all  scales.  In 
most  cases  the  vertical  velocity  spectra 
have  a  range  of  wavenumbers  with  a  -5/3 
slope,  but  the  horizontal  velocity  spectra 
do  not.  The  horizontal  spectra  appear  to 
have  ‘bumps’.  The  magnitudes  of  these 
bumps  are  emphasized  in  the  inserts  that 
have  a  linear  vertical  axis.  The  existence 
of  such  bumps  has  been  observed  in 
previous  measurements  (e.g.  refs.  19,  20, 
3)  and  explained  theoretically  in  refs.  5  & 
10.  They  are  attributed  to  a  “bottleneck” 
that  occurs  at  the  transition  between  the 
inertial  and  dissipative  range  of  the 
turbulence  spectrum. 


10  100  1000  100  1000 
k1  (rad/m) 


Figure  3.  Mean  spatial  energy  spectra.  Solid  lines: 
EJkJ,  Dashed  lines:  3/4En(k ).  Inset  figures 
are  spectra  of&^k^Eu(k^).  a  -  f  correspond  to 
Runs  A  -  F. 


Since  the  spectra  of  the  vertical  velocity  appear  to  have  a  more  extended  “inertial”  range,  we  use 
E33(kj)  to  estimate  the  dissipation  rate,  denoting  this  estimate  as  sLF.  The  results  are  summarized  in 
Table  2.  Having  data  that  extends  to  wavenumbers  in  the  dissipation  range,  enables  us  to  obtain  more 
“direct”  estimate  for  the  dissipation  (ref.  4)  using  all  or  some  of  the  measured  velocity  gradients. 
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Values  are  also  presented  in  Table  2,  along  with  the  Kolmogorov  scale,  r\  =  (v 3  /s^)14. 

Errors  associated  with  the  data  being  under  resolved  are  in  the  25%-45%  range.  There  are  significant 
spatial  variations  between  Instantaneous  distributions  comparing  eD  ,  Eqx  and  of  the  same  data 

show  but  when  averaged,  eD  and  s&  are  consistently  very  close  to  each  other  (supporting  the  validity 
of  results  that  would  be  obtained  from  vertical  profiling),  whereas  s &  is  consistently  55-64%  smaller 
than  eD  .  The  consistently  lower  values  of  seem  to  be  associated  with  the  anisotropy  of  the  oceanic 
turbulence.  The  trends  of  eLF  differ  significantly  from  the  other  methods  and  do  not  provide 
meaningful  results  for  Runs  C-F,  the  moderate  and  weak  flows. 
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One  can  characterize  the  turbulence  using  the  Taylor  micro-scale  Reynolds  number,  Re"  «  u'X/v  , 

Xu  ~  a'(\ 5v  /s)12.  To  estimate  u'  and  w'  without  being  contaminated  by  waves  one  can  use  the  second 
order  structure  function,  a  method  introduced  by  Trowbridge  (22),  and  implemented  using  PIV  data  in 
refs.  12  &  13.  Averaged  values  of  17'  and  w'  and  the  corresponding  Rex  are  also  presented  in  Table  2. 
In  Runs  A  and  B  Rex  is  in  the  300-400  range,  i.e.  a  “high”  turbulence  level.  However,  Rex  of  the 
LEO-15  site  are  68-83  at  moderate  flow,  and  14-27  when  the  mean  current  is  weak.  Since  Runs  C-F 
represent  typical  calm  weather  conditions  near  LEO-15,  the  results  indicate  that  turbulence  in  the 
coastal  bottom  boundary  layer  in  calm  weather  has  a  very  low  Reynolds  number.  The  moderate  and 
low  Re}  cases  fall  in  the  range  where  assumptions  of  universality  of  the  energy  spectrum  are  invalid 
(ref.  16). 

Table  2.  Dissipation  rates,  Kolmogorov  Scale  (r\),  vector  spacing  (8)  and  SGS  dissipation  (sSG). 
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PIV  data  can  also  be  used  for  evaluating  models  of  the  sub-grid  scale  (SGS)  stresses,  and  for 
estimating  the  SGS  dissipation  (or  energy  flux)  for  Large  Eddy  Simulations  (LES)  (Liu  et  ah,  1994).  In 
LES  the  Navier  Stokes  Equations  are  filtered  spatially  at  a  scale  A  and  the  resulting  SGS  stresses, 
xjjGS  =  Uj uj  -  UiUj  (“~”  indicates  spatial  filtering)  must  be  modelled  (see  example  in  Figure  4). 
Associated  with  this  stress  is  the  SGS  dissipation  rate,  eSG  =-xf}SSij,  Sy  =  0.5 (dujdxj  +duJ/dxj), 
which  represents  transfer  of  energy  from  the  resolved  to  subgrid  scales  Consequently,  attempts  to 
model  the  SGS  stresses  frequently  focus  on  reproducing  the  correct  levels  of  eSG .  Unlike  viscous 
dissipation,  eSG  can  be  both  positive  and  negative.  A  positive  value  indicates  flux  of  energy  from  large 
to  small  scales  whereas  a  negative  value  indicates  “backscatter”  of  energy  from  small  to  large  scales. 
Reviews  are  available  in  refs.  7,  1 1  and  14.  When  the  filter  is  in  the  inertial  range  of  isotropic, 
homogeneous  turbulence,  the  mean  eSG  is  (almost)  equal  to  the  viscous  dissipation  rate.  However,  ref. 
15  show  that  near  the  wall  of  channel  flows,  the  mean  eSG  is  small  and  even  negative. 
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Figure  4.  Instantaneous  distributions  of 
the  velocity  field  shown  in  Figure  2 
filtered  at  A -86  and  the  corresponding 
distributions  of  x13,  contoured  at 
intervals  of  5x1  O'6  nils2.  Negative  values 
are  shaded  gray. 


The  Smagorinsky,  eddy  viscosity  model  (8,  21)  for  the  SGS  stress  tensor  is  xjj  =  -2(c5A)2|)5|S)/-,  where 
|  S  1=  ^2SjjS~-  ,  A  is  the  filter  scale,  and  cs  is  a  coefficient  that  has  to  be  determined.  The  coefficient  is 
calculated  by  matching  the  measured  and  modelled  sSG ,  using  (12) 

eSG  =  ~xijSij  ~_2^'11^11  +T33^33  -t11^33  —  x 33*^11  +  12t13S13) 

We  also  estimate  the  coefficients  of  the  non-linear  model  (refs.  1,  6). 

Measured  averaged  values  of  sSG  are  presented  in  Table  2  for  three  filter  sizes,  A=4,  8  and  16  vector 

oz1  - 

spacings  (8,  Table  1).  At  high  flow  (A-B)  e  is  of  the  same  order  as  sD  ,  whereas  in  the  moderate  and 
weak  mean  flow  conditions,  the  SGS  dissipation  is  more  than  an  order  of  magnitude  lower  than  eD  , 
decreasing  to  negligible  levels  as  the  flow  diminishes.  Considering  that  cases  C-F  represent  typical 
calm  weather  conditions  near  LEO  15,  the  results  bear  significant  implications  for  applications  of  LES 
to  coastal  flows.  Ref.  12  contains  a  discussion  on  differences  between  positive  and  negative  energy 
fluxes  and  correlations  between  modeled  and  measured  SGS  stresses.  The  Smagorinsky  model 
coefficients  (Table  2)  decrease  with  decreasing  mean  flow  rate,  and  cannot  be  determined  for  run  F, 

due  to  a  negative  e  .  The  non-linear  model  coefficients  are  meaning  less  due  to  the  low  e  . 

IMPACT 

The  present  measurements  and  analysis  enable  us  to  characterize  the  flow  structure  and  turbulence  in 
the  bottom  boundary  layer  of  the  coastal  ocean  in  great  detail.  The  measurements  provide  distributions 
of  Reynolds  stresses,  spectra,  dissipation  rate,  characteristics  of  SGS  stresses  and  energy  transfer 
across  scales,  all  free  of  contamination  by  surface  waves.  These  data  are  essential  for  modeling  of 
coastal  circulation,  sediment  and  pollutant  transport,  and  biological  processes. 

TRANSITIONS 

During  FY  2000  &  2001  the  submersible  PIV  system  has  been  used  extensively  at  NSWC/Carderock 
to  measure  the  flow  structure  within  wakes  behind  several  maneuvering  models. 

RELATED  PROJECTS 

The  research  involving  measurements  of  flow  structure  and  turbulence  in  the  coastal  ocean  using  PIV 
has  been  funded  in  part  by  ONR  (present  project)  and  in  part  by  NSF  (9871961).  The  NSF  sponsored 
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project  focuses  with  the  flow/turbulence  in  the  water  column  above  the  boundary  layer,  whereas  the 
present  project  focuses  on  the  boundary  layer. 
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